Electrostatic fluctuations promote the dynamical transition in proteins 
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Atomic displacements of hydrated proteins are dominated by phonon vibrations at low temper- 
atures and by dissipative large-amplitude motions at high temperatures. A crossover between the 
two regimes is known as a dynamical transition. Recent experiments indicate a connection between 
the dynamical transition and the dielectric response of the hydrated protein. We analyze two mech- 
anisms of the coupling between the protein atomic motions and the protein-water interface. The 
first mechanism considers viscoelastic changes in the global shape of the protein plasticized by its 
coupling to the hydration shell. The second mechanism involves modulations of the motions of 
partial charges inside the protein by electrostatic fluctuations. The model is used to analyze mean 
square displacements of iron of metmyoglobin reported by Mossbauer spectroscopy. We show that 
high flexibility of heme iron at physiological temperatures is dominated by electrostatic fluctuations. 
Two onsets, one arising from the viscoelastic response and the second from electrostatic fluctuations, 
are seen in the temperature dependence of the mean square displacements when the corresponding 
relaxation times enter the instrumental resolution window. 
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I. INTRODUCTION 

Measurements of the Mossbauer absorption of 57 Fe in 
metmyoglobin crystals revealed that the mean square dis- 
placement (msd) of this atom starts to grow faster with 
increasing temperature above Try ~ 200 This find- 
ing was followed by similar observations from neutron 
scattering,— which by now have been reported for a large 
number of proteins and other biopolymers,— all demon- 
strating the same phenomenology^^ The increase in the 
slope of the protein msd as a function of temperature was 
called a "dynamical transition" , presently assigned to a 
rather broad range of onset temperatures Try — 200 — 240 
K. The basic observation is that the high-temperature 
flexibility of the proteins much exceeds the linear extrap- 
olation of the low-temperature behavior characteristic of 
an expanding solid. The low-temperature msd is well 
characterized by the observed phonon spectrum of the 
protein,— while the high-temperature msd excess is linked 
to dissipative long- wavelength modes with energies below 
~ 3 meV^iS 

Early explanations of the dynamical transition offered 
scenarios ranging from detrapping of the protein con- 
formational motions from low-energy states^ ' 11 ' 12 to the 
glass transition in bulk water^ Several recent observa- 
tions have shifted the focus to the protein hydration 
shell. The disappearance of the dynamical transition in 
dry protein powders^ and its separate existence for the 
hydration water ; 14 ' 15 point to a strong link between the 
dynamical transition and the protein hydration shell. 

Despite somewhat different semantics, the views in the 
field seem to converge to the notion of the critical role 
of the hydration shell in driving the dynamical transi- 
tion. According to Doster^ "The onset of the dynamical 
transition depends on the solvent viscosity near the pro- 



tein surface.. . .The protein- water a-process consists of a 
concerted librational motion of protein surface residues, 
coupled to translational jumps of water molecules on the 
same time scale." The physical mechanism behind the 
transition is assigned in this view to the caging of the 
protein by water's hydrogen bonds, stiffening its confor- 
mational flexibility. As the temperature increases, the 
population of broken hydrogen bonds grows exponen- 
tially, resulting in a release of the protein conformational 
flexibility in a narrow range of temperatures. Although 
appealing, this concept does not address the key question 
of how ~ 0.5 — 2 ps events of hydrogen bond breaking 
develop into a ~ 2 ^ts collective relaxation process at Try. 

Fraucnfclder and co-workers have recently suggested a 
somewhat different scenario, also focusing on the dynam- 
ics of the protein hydration shclli 1 ^— According to their 
view, protein conformations are coupled to motions of 
the hydration layer, with a relaxation time following an 
Arrhenius law. This relaxation mode is therefore con- 
sidered to be a secondary, /3-process according to the 
common classification adopted in the field of supercooled 
liquids 3 This secondary relaxation is assessed by dielec- 
tric spectroscopy of protein samples embedded into solid 
poly(vinyl)alcohol (PVA), thus eliminating the bulk wa- 
ter relaxation from the dielectric response! 18 ' 19 

The use of dielectric absorption of PVA-confincd pro- 
teins yields a surprisingly accurate account of the temper- 
ature dependence of the Lamb-Mossbauer factor equal to 
the fraction of recoillcss absorption, / = exp[— kg((Sx) 2 )]. 
Here, ((Sx) 2 ) is the msd of the heme iron of metmyo- 
globin in projection on the wavevector ko of 7-radiation. 
The iron msd can be separated into a vibrational, low- 
temperature component ((Sq) 2 ) oc T, described by the 
vibrational density of states (VDOS) of the protein, 
and a high-temperature component ((SQ) 2 ) appearing at 



2 



T > Tjj. Correspondingly, / = / 9 /q becomes a product 
of two components, f q and fq. It turns out that the 
temperature variation of the Lamb-Mossbauer factor fq 
originating from the high-temperature msd is exception- 
ally well described by the variance of the sample dipole 
moment at the same hydration level 

f Q = cM-k 2 {(SQ) 2 )} = ((5Mf) < /((6Mf). (1) 

In this equation, ((<5A/) 2 )< is defined as the integral of 
the frequency-dependent variance of the sample dipole 
moment M u over the frequencies below (subscript "<") 
the instrumental frequency w bs = l/^obs, ^obs = 140 
nsJ£ The parameter fq then determines the fraction of 
the sample dipole that has not had a chance to alter on 
the life-time of the iron nucleus, thus keeping the nuclei 
in resonance for Mossbauer absorption. 

The empirical connection between the dynamical tran- 
sition and dipolar fluctuations is additionally supported 
by recent observations of breaks in the dependence of the 
terahertz dielectric absorption on temperature at typi- 
cal values of Tdj21 All these observations, although ad- 
vancing the field toward identifying the physical modes 
responsible for high-temperature flexibility of proteins, 
pose a significant conceptual challenge. 

Both Mossbauer and neutron-scattering techniques 
probe translational atomic motions on their correspond- 
ing resolution windows. It seems therefore natural to 
relate the break in the temperature dependence of the 
msd to changes in the dynamic and/or static properties 
of atomic translations. This is the conceptual frame- 
work behind the glass-transition scenario^ which, even 
in the current form emphasizing the hydration layer, is 
focused on the caging arrest, i.e., on the primary effect 
of short-ranged repulsive interactions in the system. On 
the contrary, the dielectric measurement a 19 ' 21 i 22 shift the 
focus to the long-ranged electrostatic interactions, which 
is what dielectric spectroscopy is sensitive to at the first 
place. What is the correct view? 

Answering this question requires gaining deeper in- 
sights into the actual physical modes coupled to the pro- 
tein msd at high temperatures, the problem that has 
evaded direct experimental inquiry so far. Numerical 
simulations point to enhanced fluctuations of hydrogen 
bonds of hydration water at high temperatures^ but 
those can be projected on either density or dipolar col- 
lective modes. Since a strong link between the dynam- 
ical transition and the hydration shell has been clearly 
established, the main question posed by recent studies 
is whether density or orientational collective modes drive 
the transition^ They are mostly decoupled by symmetry 
and can therefore be considered as two distinct mecha- 
nisms of altering the protein flexibility. 

The goal of this paper is to present some initial esti- 
mates of the relative importance of the density and ori- 
entational fluctuations in the protein dynamical transi- 
tion. We model the density fluctuations at the protein- 
water interface by viscoelastic response and the orien- 
tational dipolar fluctuations in terms of electrostatic re- 



sponse. The dependence of the onset temperature Tr> on 
the observation window is an important ingredient of the 
observations ; 16 i 25 which is introduced into the theory by 
limiting the range of frequencies over which the response 
functions are integrated^ similarly to Eq. ([1]). We start 
with formulating the model, followed by the results of 
calculations. 



II. MODEL 

The purpose of our model is to determine the msd of 
a single atom, heme iron of metmyoglobin, as a func- 
tion of temperature. The internal motions of the protein 
can roughly be separated into two modes, the phonon 
vibrations q and dissipative large-scale motions Q. Even 
though phonons do not formally exist in proteins, we will 
use this language to distinguish short-wavelength vibra- 
tions from motions altering the protein's global shape. 
Accordingly, we will split the coordinate of the iron atom 
r = q+Q into two statistically independent components, 
q and Q. The former can be expanded in protein's nor- 
mal modes or directly calculated from the VDOS mea- 
sured, for instance, by the phonon-assisted Mossbauer 
scatteringi 6 ' 27 i 28 The corresponding msd of the vibra- 
tional coordinate q is, in the classical limit, a linear func- 
tion of temperature 

<(^) 2 } = a q T. (2) 

The proportionality coefficient a q is calculated from the 
VDOS according to the standard prescriptions. 

The dissipative motions of the protein are seen in neu- 
tron scattering spectra as a quasielastic peak with ener- 
gies below ~ 4 meV, growing in intensity with increasing 
temperature £ With the sound velocity in a protein^^ 
of ~ 1700 m/s, the wavelength of the corresponding vi- 
brations is about 26 A, which is comparable with the 
diameter of myoglobin, 2R = 36 A. These modes there- 
fore alter the global shape of the protein, which is the 
domain of the viscoelastic response. 

In order to obtain a first-order estimate of the geom- 
etry changes involved, we will model the protein mo- 
tions as radial viscoelastic vibrations of a sphere of ra- 
dius R immersed in a viscoelastic water continuum. The 
low-frequency variance of the iron coordinate is then 
simply related to the radius fluctuations of the sphere 
as ((<5Q 2 )} = (r/R) 2 ({5R) 2 ). The latter can be found 
by solving the standard equations of viscoelasticit y 30 ' 31 
yielding the response function Xfl( w ) connecting the 
change of the sphere's radius R to an oscillatory pres- 
sure p(t) = p e lut applied to the sphere's surface. The 
result is2£ 

Xi?(w) = ~4irR3AK p (u)+^ w (u)' (3) 

Here, AK p (u>) = K p {tu) — K P} q is the viscoelastic bulk 
modulus of the protein minus its bulk modulus K Pi q a t 
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zero frequency. Further, fi w {ui) is the shear modulus of 
water. Applying the fluctuation-dissipation theorem^ 2 - to 
Eq. one gets 

((SQ.) 2 ) = -g^Im 1 (4) 

where V p is the protein volume. 

We now proceed to calculating the Lamb-Mossbauer 
factor^! 



/=(|(e*>*>| 2 > het , 



(5) 



where x = ko • (q + Q) is the projection of the iron 
displacement on the direction of photon propagation, 
k = ko/fco- There are two averages in this definition: 
the inner angular brackets denote an ensemble average 
over the protein and water modes affecting the position of 
iron in a single protein, while the outer angular brackets 
denote the average over the proteins in the sample. This 
second average carries the subscript "het" to emphasize 
that it reflects the heterogeneity of the sample, such as 
for instance variations in the hydration level among dif- 
ferent proteins in the protein powder. We do not consider 
the heterogeneous average in our present study and limit 
ourselves by the inner average only. This approxima- 
tion amounts, in experimental techniques, to considering 
the narrow feature of the Mossbauer absorption line and 
subtracting the broad base-line originating from the sam- 
ple heterogeneity— Accordingly, the experimental points 
shown in Fig. [1] are obtained from the area f(T) of the 
narrow line as - In f(T)/(k ) 2 , fc = 53.2 A -1 . 

The ensemble average over the water/protein statis- 
tical distribution can be described in terms of the free 
energy F(x) such that the inner brackets in Eq. ([5|) be- 
come 



(e lk ° x ) 
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where & = l/(k B T) and Z = J cx-p{-f3F(x)]dx. The free 
energy F(x) is determined by projecting^ the manifold 
of q and Q coordinates on the single coordinate x 



-f3F(x) 



S[x - k • (q + Q)] 

/ p -/3ffo(q,Q)-/3z<Mq,Q)\ 



(7) 



dqdQ. 



In this equation, Ho(q, Q) is the Hamiltonian of classical 
harmonic vibrations of the protein and <j) w is the elec- 
trostatic potential of the surrounding dielectric medium 
acting on the heme iron carrying charge z. The subscript 
"w" in Eq. ([7]) specifics water as the main source of the 
electrostatic fluctuations. It is not, however, required by 
the theory, and slow protein motions, not included in the 
calculation of ((5q) 2 ), can contribute to the fluctuations 
of the electrostatic potential as well (see below). 

The Hamiltonian Ho(q, Q) can be given in the Gaus- 
sian form 



5q 2 



SQ 2 



2((Sq) 2 ) 2((SQ) 2 ) : 



(8) 



where the variance of q is given by Eq. <j2j) and the vari- 
ance of Q requires additional explanation. 

The limited instrumental time r b s affects the observ- 
ables and, in fact, the dynamical transition itself becomes 
possible only when the relaxation time of a collective 
mode of the hydration shell coupled to high-temperature 
protein's motions enters the experimental observation 
windowi 16 ' 19 Therefore, the variance of the slow disper- 
sive motions of the heme iron is not a thermodynamic 
variable referring to an infinite observation window, but 
a property affected by instrumental resolution! 26 i 36 This 
is reflected by the subscript ">" in Eq. ([5]) which indi- 
cates that ((SQ) 2 ) > is calculated by integrating the re- 
sponse function in Eq. (f4|) over the frequencies exceeding 
the observation frequency w b s 
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((SQ) 2 )> 



({6Q„) 2 ){du/w). 



(9) 



The statistical average over the electrostatic fluctua- 
tions can be simplified by a first-order expansion of the 
potential 4> w in x: </)«,(q, Q) — 4> w ,o ~ xE W) where 4>wfl 
is the potential at the equilibrium position of iron and 
E w is the electric field projected on ko. Assuming that 
E w is a Gaussian variable, one gets a Gaussian form of 
[3F(x) = x 2 / {2{{5x) 2 )) where the variance {(Sx) 2 ) is 



((Sx) 2 ) = ((Sx) 2 ) el /M E 



(10) 



Here, the elastic msd ((6x) 2 ) cl = ({Sq) 2 ) + ((<5Q) 2 }> is 
the sum of two statistically decoupled components. Fur- 
ther, the correction Me represents softening of atomic 
vibrations by electrostatic fluctuations of the hydration 
shell. 

Similarly to the viscoelastic effect, the electrostatic 
softening depends on the observation window. Account- 
ing again for the frequency cutoff introduced by the finite 
instrumental resolution, it is given in the form 

/>oo 

M E = l-(f3z) 2 ((6x) 2 ) cl C E (uj)(kj/(2n), (11) 

where Ce(u;) is the Fourier transform of the time auto- 
correlation function of the field E w (t) 

{5E w {t)5E w {0))e lult dt. (12) 

-oo 

By applying the fluctuation-dissipation theorem^ once 
again, one can recast Eq. (|11| in terms of the response 
function xe{u) representing the polar response to an os- 
cillating dipolc probe m(t) = mo exp(iojt) placed at the 
position of the iron atom 



M E = 1 - I3z 2 ((6x) 2 ) cl I XB(w)dw/(7rw). (13) 



OO 



This form of the correction factor accounting for the dipo- 
lar fluctuations of the water shell is used in the calcula- 
tions below. 
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FIG. 1. Temperature dependence of the msd of iron in met- 
myoglobin. Experimental msd^ are shown by circles. Di- 
amonds refer to the msd calculated from the VDOS mea- 
sured from phonon-assisted Mossbauer effect at different 
temperatures,— whereas the dashed line is the same calcula- 
tion from the VDOS at 235 KM The dash-dotted line refers 
to the msd combining phonons with viscoelastic oscillations of 
the protein shape and the solid line refers to the combination 
of all three effects: protein's phonons, viscoelastic shape oscil- 
lations, and dipolar fluctuations of the hydration layer. The 
solid line is obtained with AK P (T) as explained in the text, 
while the dotted lines refer to the temperature-independent 
AK P = 3.2 GPa obtained from the fit of Eq. (Tl7j) to the 
protein intrinsic compressibility at 298 K.— The inset shows 
the rise of the msd due to the elastic motions near the glass 
transition temperature of the protein, T g ~ 180 ± 15 K. 



III. MEAN SQUARE DISPLACEMENT 

Here we outline the calculations performed using the 
model developed in this paper. We need to mention 
that many parameters entering the model are not ex- 
perimentally available. Some of them can be potentially 
extracted from numerical simulations. The usefulness of 
simulations is, however, limited for interpreting the ex- 
perimental data since reproducing heterogeneous condi- 
tions of partially hydrated protein powders presents sig- 
nificant challenges to simulation protocols. Likewise, the 
viscoelastic model used here should be viewed as only a 
first step toward a more realistic description of the clastic 
response of hydrated proteins. However, one of the major 
conclusions of this paper is the dominance of electrostat- 
ics in the high-temperature flexibility of proteins and a 
relatively small effect of viscoelastic motions on the iron 
msd. This observation puts high priority to the develop- 
ment of the electrostatic component of the model, and 
makes the limitations of the modeling of the viscoelastic 
response less critical. 

The viscoelastic response functions entering Eq. ([3]) 
were taken in the Maxwell forrn^ 



AK p iuiTp(T) 
1 + iuT p (T) 
G 

CO 

iur w (T) 



(14) 



iuir w (T) ' 



In this equation, AK p — K poc — K p _q is the change in 
the bulk protein modulus between infinite and zero fre- 



quencies and Goo is the high-frequency shear modulus 
of water; t P)W (T) are the corresponding relaxation times. 
From the Maxwell equation, one gets the shear viscosity 
Vw(T) = Goo(T)t w (T) which is well tabulated down to 
the water nucleation temperature^ The shear relaxation 
time was obtained from r] w (T) and Goo(T)/GPa = 1.68 — 
0.0127(T - 273) taken from Ref. [H The resulting t w (T) 
turns out to be close to the exponential relaxation time of 
the longitudinal modulus extracted from inelastic x-ray 
scattering^ n(T)/s = 0.84 x lCT 15 exp(1910 K/T). 

The VDOS of metmyoglobin powders is well estab- 
lished by phonon-assisted Mossbauer measurements fi&iSL 
and the temperature slope a q in Eq. ([2]) was calculated 
from the experimental VDOS (dashed line in in Fig. [IJ. 
Figure Q] also shows a rise in the msd due to the onset of 
viscoelastic oscillations of the protein (dash-dotted line) . 
These results were obtained by adopting Eq. (fT4|) for the 
elastic moduli with the protein's relaxation time t p (T) 
from the measurements done on dry protein powders^ 
The relaxation process in dry proteins is too slow to 
enter the observation window of the spectrometer and 
AK p (uj) ~ AK p in Eq. (fT4"j) . 

This notion implies that that the frequency depen- 
dence of the moduli, e.g., Debye vs dispersive relaxation 
does not significantly affect the outcome of the calcula- 
tions and only AK p (T) matters for the temperature de- 
pendence of the viscoelastic msd. The latter was adopted 
to reproduce the experimental intrinsic compressibility 4 ^ 
of myoglobin 3 ! j3 T = 11.04 Mbar" 1 at T = 298 K and the 
temperature variation of Young's moduli of myoglobin 
crystals at lower temperature o 42 ' 43 (see below). The use 
of experimental compressibility to parametrize the model 
also implies that the overall viscoelastic component of 
atomic displacements is limited by the thermodynamic 
experimental value and can only be lower for a given in- 
strumental observation window. 

The most uncertain part of our analysis is the calcu- 
lation of the electric field response function xe(w). It 
can be calculated by solving the Poisson equation for the 
heme immersed in the heterogeneous dielectric formed 
by the protein and its hydration layer. However, the as- 
signment of the dielectric constants to both the protein 
and the thin layer of water surrounding it in the powder 
is subject to significant uncertainties. We will therefore 
restrict ourselves to rough estimates aimed to establish 
whether electrostatic fluctuations can produce a signifi- 
cant effect on the msd under a reasonable set of assump- 
tions. The results of fitting are additionally supported 
by Molecular Dynamics simulations of the fully hydrated 
metmyoglobin^ 

The dielectric response to charges immersed in a po- 
lar medium is dominated by longitudinal modes of dipo- 
lar polarization 4 ^ characterized by the dielectric modulus 
e(u) , where e(w) is the complex, frequency-dependent 
dielectric constant of the protein-water mixture. The 
dielectric response function \e{^>) establishes the reac- 
tion field of the dielectric medium in response to a probe 
dipole placed at the position of the heme iron. It scales as 
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FIG. 2. Dielectric longitudinal relaxation time t^(T) = 
(^cx,/^ s )te(T) reported in the literatur e) 19 ' 47 ' 48 with the pro- 
tein hydration levels indicated in the legend. The results from 
Ref.— refer to myoglobin confined in PVA, results from Ref. 24 
are for metmyoglobin crystals, and Refs . 47 ' 48 refer to met- 
myoglobin powders. The solid line refers to the longitudinal, 
Debye relaxation time of the reaction field response function 
X e(u) [Eq. |[I5)]. t l {T)/s = 3.2 x 1CT 13 exp[3000 K/T], used 
to fit the experimental msd for the heme iron.— 



the inverse cube of the characteristic size d of the heme. 
The loss function Xe^) can therefore be written in the 



form 



Xe[uj) - & WW 



(15) 



Dielectric properties of partially hydratcd proteins 
have not been well characterized since the results are 
strongly affected by both the sample preparation and 
the hydration level. Even dielectric relaxation times mea- 
sured on samples of close hydration level are rather incon- 
sistent. This point is illustrated in Fig. [2] where results on 
partially hydratcd myoglobin powders of hydration level 
h — 0.3 — 0.5 (in g of water per g of protein) have been 
assembled ! 19 ' 46 " — We also show measurements done on 
myoglobin crystals 2 ^ for the sake of comparison. Multiple 
relaxation processes are common for such measurements, 
and the fastest relaxation, commonly attributed to the 
hydration shell , 22 ' 47 is shown in Fig. El 

Even if we could firmly establish the proper relaxation 
time for the the sample dipolc moment, this would not 
necessarily give us the relaxation time of the reaction 
field correlation function required for x'ki^)- Because 
of the linear scaling of the dipole moment variance with 
the number of dipoles, dielectric measurements empha- 
size the effect of outer solvation shells, while Xb( w ) 1S 
dominated by waters closest to the probe dipole (heme's 
iron). In view of these uncertainties, we have constructed 
the temperature-dependent relaxation time te(T) that, 
together with the parameter d in Eq. (fT5j). allows us to 
fit the experimental msd. 

Assuming the Debye form for the fast relaxation com- 
ponent in e(a;) j 22 ' 48 the response function in Eqs. (fj"3)) 
and (|15|) gains the form 
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<5 2 i + {uT L y 
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FIG. 3. Intrinsic isothermal compressibility Pt{T) of a pro- 
tein calculated from Eq. (|17[) with uj D hs = 1 (dashed line), 
10 2 (solid line), and 10 3 MHz (dash-dotted line). The protein 
and water parameters are those adopted for the calculation of 
the myoglobin msd, with AK p (T)/GPa, = 3.22 - 0.03 X (T — 
298) + 0.00025 x (T - 298) 2 obtained to fit the protein intrinsic 
compressibility at 298 K 3 ^ and the temperature dependence 
of Young's moduli (crosses in the inset). The inset shows 
experimental expansivity of the hydration shell of lysozyme 
(circles) j£i experimental inverse Young's moduli of myoglobin 
crystals exposed to air of 95-100% (diamonds)^ 3 - and 75% hu- 
midity (crosses).— The solid line shows Pt{T) calculated from 
Eq. (| lTp at u bs = 10 2 MHz; all curves are normalized to the 
corresponding values at T = 298 K. 



where the parameter 5, 5 
tecteristic length and Cq = 



(3z 2 co/d 3 sets up a chara- 
L — ej 1 is the Pekar factor. 
In addition, t£ = (foo/ '^s)te is the longitudinal dielectric 
time and eoo and e s are the high-frequency and static di- 
electric constants, respectively. The parameters 5 = 0.12 
A and t l (T)/s = 3.2 x 10~ 13 exp[3000 K/T] (solid line 
in Fig. [2|) were used in fitting the experimental msd. 
The fitting relaxation time tl(T) is generally consistent 
with dielectric measurements and, in addition, the Arrhc- 
nius slope of t l (T) matches our simulations of the pro- 
tein Stokes-shift dynamics at elevated temperatures^ 
More detailed calculations might require replacing one- 
relaxation Debye dynamics in Eq. (|16p with dispersive 
dynamics characterized by a distribution of relaxation 
time, as suggested by the NMR experiment.— 

The overall fluctuations of the protein volume are de- 
termined by the response function in Eq. (|3|) combining 
the dynamic elastic moduli of the protein and the hydra- 
tion shell. This connection can be used to parameterize 
the model on volumetric properties of hydrated proteins, 
in particular on protein's intrinsic compressibility4i For 
a given instrumental resolution, one obtains from Eq. ([3]) 
for the isothermal compressibility /3t °c ((3Vp) 2 ) of the 
protein 

/>oo 

(i T = -(6/tt) / Im [3AA p (w) + W^f 1 (du/u). 



(17) 

In Fig. [3] we show (3t(T) for the parameters adopted 
in the calculations of the iron msd and several values of 
u; b s - The intrinsic compressibility of the protein rises 
sharply at the point close to protein's glass transition 
T g . The latter depends on the observation window, but 
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FIG. 4. The loss function /32 2 Xb(^) obtained from the fitting 
of the experimental msd to Eqs. (|10p . (| 13p . and (|16|l (marked 
as "F") and from direct MD simulations of the electric field 
acting on the iron of metmyoglobin. The plot shows the result 
for the overall electric field produced by protein and water 
("P+W") and by protein ("P") and water ("W") separately. 
The simulation trajectories (T = 300 K) were 45 ns long, 35 ns 
of which were used for collecting the correlation functions—^ 

is close to reported values T g ~ 180 ±15 K 43 i 51 i 52 marked 
by breaks in several observable parameters^* 2 - The rise 
of compressibility at T g in our calculations is caused by 
the water component of the viscoelastic response func- 
tion when the relaxation time t w becomes smaller than 
T b s . This result is consistent with the glass transition 
of the hydration shell expansivity^ shown in the inset in 
Fig. [3] The inset in Fig. [3] also shows compressibilities 
obtained from experimentally reported Young's moduli 
of myoglobin crystal s 42 ' 43 assuming that their Poisson 
ratios are independent of temperature. The temperature 
variation of AK P (T) in our calculations shown in Fig. [1] 
was chosen to match these data. 

The fit of AAp(T) to crystalline Young's moduli 
results in an upward increase of the elastic msd at 
the highest temperatures shown in Fig. [TJ This up- 
ward increase reflects pre-melting of myoglobin crys- 
tals when their Young's moduli approach zero.— Since 
the melting temperature is typically higher in protein 
powders,— AK p (T) obtained from fitting the crystal 
data might overestimate these effects; the iron msd with 
a temperature-independence A.K p = 3.2 GPa is shown 
by the dotted line in Fig. [TJ 

The results shown in Fig. [TJ indicate that electrostatic 
fluctuations far outweigh viscoelastic vibrations in the 
iron msd. We additionally confirm this outcome by com- 
paring the function /3z 2 x^(w) from our fitting to the 
same function obtained from MD simulation of the fully 
hydratcd metmyoglobin.— Figure [4] shows the response 
functions from the electric field fluctuations produced by 
the protein and water combined and by each component 
separately. The height of the maximum quantifies the 
strength of the msd modulation by the corresponding 
electrostatic component, and it is of main importance 
for this comparison. 

We have assumed so far that the protein is electrostati- 



cally non-polar, and its hydration shell is the main source 
of the electrostatic fluctuations. It does not need to be 
so. Low-frequency motions, not included in the VDOS 
used to caculate ((5q) 2 }, can modulate the protein's par- 
tial charges (dipole moments of a-helices, ionized sur- 
face residues, etc.) and compete in the elecrostatic noise 
with the hydration layer. This might be particularly 
true for partially hydratcd protein powders where the 
fluctuations of the water dipoles are probably reduced 
to motions of polarized domains around ionized surface 
residues. 

Figure [U in fact shows that the protein component of 
~x!e(oj) exceeds that of water, and its maximum is higher 
than that of the fitting function from Eq. (fT6|) . The elec- 
trostatic fluctuations of the protein itself are therefore 
sufficient to produce the observable msd and, in addition, 
our estimates do not seem to overestimate the effect of 
the electrostatic fluctuations on the msd. The primary 
role of water in powders might be reduced to ionizing 
the surface residues of the protein and plasticizing its 
motions above T g (Fig. [3|). Water in patches solvating 
ionized residues is strongly coupled to the protein both 
electrostatically and by surface hydrogen bonds. The re- 
laxation times of their electrostatic response functions 
are therefore close (Fig. 3]), resulting in matching on- 
set temperatures of the dynamical transition for each 
component i 14 ' 15 

Further, the overall loss function Xjj(w), which in- 
cludes cross-correlations between the water and protein 
electric fields, shows a slower relaxation time than each 
component separately. The relaxation time of 6.3 ns of 
the essentially Debye overall function Xb( w ) i s dose to 
tx(300 K) ~ 7 ns adopted in fitting of the experimental 
msd. It is this loss function, combining the protein and 
water electrostatics, that is of primary interest for the 
modeling of the high-temperature flexibility of proteins. 

IV. DISCUSSION 

The picture presented here assigns an increase in the 
protein msd at the dynamical transition to the en- 
trance of a collective relaxation time of the protein- 
water interface into the observation window of the 
spectrometer j 19 i 22 i 53 We consider two types of interfacial 
fluctuations, elastic modes changing the global shape of 
the protein and electrostatic fluctuations. Electrostat- 
ics turn out to be the main factor affecting the high- 
temperature portion of the msd. 

The longitudinal relaxation time of the electric field 
fluctuations, tl(T), determines the transition tempera- 
ture by the condition lu i, s tl(Te)) ~ 1. With the Arrhe- 
nius form for the relaxation time tl(T), this condition 
predicts a logarithmic dependence of Tjj on the observa- 
tion frequency, 

T D cx lln^obsTojr 1 , (18) 
where tq is the preexponent in tl(T). For instance, with 
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FIG. 5. Variance of the water's electrostatic potential at the 
active site, Aq 2 ((S(f>) 2 ) of the protein plastocyanin from MD 
simulations (circles).— Diamonds show the difference of water 
potentials in equilibrium with the active site carrying charges 
qi and q%, Aq((<f>)i — {4>)2), Aq = q- 2 — q\. This latter 
quantity is sensitive to high-frequency ballistic modes of the 
hydration water, but not to collective fluctuations of the shell 
dipole.— The two calculations coincide in the linear response 
approximation, which is valid at low temperatures. Linear 
response breaks down when the collective mode of water's 
dipolar polarization enters the observation window fixed by 
the length of the simulation trajectory. The spike at ~ 220 
K in the potential variance carries signatures of a weak first- 
order transition, but its origin is currently unclear. 



the observation window of neutron scattering of ~ 500 
ps and of Mossbauer spectroscopy of 140 ns, the above 
equation yields 1.4 for the ratio of Td values measured by 
neutron and Mossbauer techniques (r = 10™ 13 s). This 
estimate assumes equal electrostatic relaxation times for 
(mostly surface) protons and heme iron, which is likely 
not true. The actual picture is also more complex as 
several slope changes contribute to the overall tempera- 
ture dependence of the msd^ It is also the case with the 
present model producing two different onsets arising from 
viscoelastic and electrostatic fluctuations. An increase in 
Td was also reported for proteins solvated in glycerol 
and in concentrated sucrose- water solutions^ Although 
an increase in viscosity does shift Td in the right direc- 
tion according to Eq. (fl8|) , the alteration of the effective 
polarity of the hydration layer and the surface charge 
distribution of the protein might be other factors con- 
tributing to the shift. Generally, the the present model 
predicts a decrease in the protein atomic displacements 
for hydration in solvents of lower polarity. 

The main physical question looming behind the phe- 
nomenon of the dynamical transition is what are the 
mechanisms and physical modes allowing high flexibil- 
ity of proteins at physiological temperatures. We em- 
phasize here electrostatic fluctuations as the primary ori- 
gin of the increase in the protein's atomic displacements. 
This mechanism connects the translational manifold of 
the protein's interior to the dipolar oricntational mani- 
fold of the hydration layer. While this connection was es- 
tablished empirically by experiment [Eq. (pQ)], numerical 



simulations directly show the same basic phenomenology 
for the electrostatic fluctuations and the atomic msd. 

Figure [5] shows the results of numerical simulations for 
the variance of the electrostatic potential produced by 
the water hydration shell at the active site of the protein 
plastocyanhij 5 - 5 - A break in the temperature dependence 
at Td refers to the time-scale of ~ 10 ns fixed by the 
length of the simulation trajectory. The difference of the 
first moments of the potential in the two redox states 
of the protein (diamonds in Fig. [5]) gives the component 
of the same property produced by the ballistic dynam- 
ics of the hydration shell and not sensitive to its collec- 
tive relaxation^ In a sense, the diamonds in Fig. [5] are 
analogs to the diamonds in Fig. [T] referring to the vibra- 
tional component of the msd£ There is a clear qualitative 
similarity between laboratory and numerical results pre- 
sented in Figs. Q] and [5] 

Because the response function of the water's electric 
field scales as d -3 with the distance d from the sur- 
face inside the protein, dipolar fluctuations of the hy- 
dration shell will mostly affect protein's surface residues. 
The vibrations of the surface protons will therefore be 
softer than of interior protons, and they will stronger 
contribute to the observable msd. There is also a pos- 
sibility of "surface melting" when Me = in Eq. (fTU|) 
is reached with rising temperature for a group of atoms. 
The low-temperature conformation of the corresponding 
residues will become unstable, with the instability re- 
leased through a conformational transition. 

The present model predicts higher flexibility for atoms 
carrying higher partial charges. In case of heme iron 
this implies higher flexibility of the protein oxidized state 
compared to the reduced state. While this prediction 
qualitatively agrees with experiment , 57 ' 58 more detailed 
studies are required to distinguish the effect of electro- 
static fluctuations from the alteration of the VDOS also 
occurring upon changing the redox state. 

V. CONCLUSIONS 

The model proposed here treats high-temperature 
atomic displacements of the protein as a combination of 
viscoelastic deformation of the global protein shape and 
electrostatic fluctuations coupled to the atomic charge. 
We suggest that electrostatic fluctuations dominate the 
high-temperature flexibility of proteins. 
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